Report for Model 9

— for both jointly and separately

2024-10-08

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      1.013       1.08
## LDevsum      1.016       1.02
## dh0          0.992       1.00
## dh1          0.994       1.01
## dl0          2.242       4.93
## dl1          1.018       1.04
## dl2          1.119       1.46
## dl3          3.968       9.00
## 
## Multivariate psrf
## 
## 2.72

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1232.80062 22803.7627
DIC3 1162.49887 20640.9861
PWAIC 43.86491 993.4114
WAIC 1190.09517 21023.4924

2. Separate Modeling of High-Level

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.01       1.03
## dh0           1.02       1.07
## dh1           1.02       1.06
## dl0           1.36       2.37
## dl1           1.80       3.49
## dl2           1.01       1.07
## dl3           2.46       4.91
## 
## Multivariate psrf
## 
## 1.89

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

H0a
DIC 1220.58844
DIC3 1158.83229
PWAIC 43.34125
WAIC 1186.21081

3. Separate Modeling for Low-level

3.1 Trace plots for convergence check

The current MCMC setting is:

  • 300000 iteration;
  • 200000 burn-in;
  • 100 thinning.

3.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## Devsum        1.01       1.05
## LDevsum       1.01       1.03
## dl0           1.37       2.49
## dl1           2.57       5.95
## dl2           1.02       1.05
## dl3           2.43       6.74
## 
## Multivariate psrf
## 
## 1.87

3.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

3.4 WAIC results

L7
DIC 23921.990
DIC3 21168.527
PWAIC 1536.806
WAIC 21915.208